seem (version 1.0)

Temp response functions: Functions to define temperature response.

Description

Functions parab, kt.rate, Q10 are for single temperature response. Function Q10.mruns explores sensitivity to parameters of Q10.

Usage

parab(Temp,a,b)
kT.rate(Temp,k0,q10)
Q10(Temp, param)
Q10.mruns(Temp, param.nom, Topt.sens, Tmax.sens, q10.sens)

Arguments

Temp

temperature

a

parameter

b

parameter

k0

parameter

q10

parameter

param

parameter array c(Topt,Tmax,q10) where Topt is temp for which optimum occurs, Tmax is maximum temperature, and q10 is coefficient

param.nom

nominal values of parameters c(Topt,Tmax,q10)

Topt.sens

values to vary parameter Topt

Tmax.sens

values to vary parameter Tmax

q10.sens

values to vary parameter q10

Value

• Temperature response

Details

The three functions parab, kT.rate, and Q10, calculate the response to temperature using the parabolic, kT by q10, and the gamma-like q10.

Function Q10.mruns explores sensitivity to parameters of Q10.

References

Acevedo M.F. 2012. Simulation of Ecological and Environmental Models. CRC Press.

See Also

Simulation functions sim, simd

Examples

# parabolic temp response
Temp <- seq(0,50); nT <- length(Temp)
Tu <- c(40,45,50);np <- length(Tu)
Tl <- 0
# dimension array
y <- matrix(nrow=nT,ncol=np)
# loop for the various optimal temperatures and then for temp range
for(j in 1:np){
 for(i in 1:nT) y[i,j] <- parab(Temp[i],Tl,Tu[j])
}

# rate kT
Temp <- seq(20,30); nT <- length(Temp)
k0 <- 1
q10 <- c(1,1.5,2); np <- length(q10)
kT <- matrix(nrow=nT,ncol=np)
for(i in 1:np){
kT[,i] <- kT.rate(Temp,k0,q10[i])/kT.rate(Temp[1],k0,q10[i])
}

# Q10
Temp <- seq(0,50); nT <- length(Temp)
Topt= 25; Tmax <- 40; q10 <- 1.5
fT <- Q10(Temp, param=c(Topt,Tmax,q10))

Temp <- seq(0,50); nT <- length(Temp)
# param Topt, Tmax,q10
param.nom <- c(25,40,1.5)
Topt.sens <- c(20,25,30)
Tmax.sens <- c(40,45,50)
q10.sens <- c(1,1.5,2)
np <- length(q10)
test <- Q10.mruns(Temp,param.nom,Topt.sens,Tmax.sens,q10.sens)